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ABSTRACT 

We report the detection of radio emission coincident with the ultraluminous X-ray source (ULX) 
in Holmberg II. The radio emission is diffuse and resolved, covering an area ~ 60 x 40 pc in extent 
and well-matched to the recently discovered He n nebula surrounding the X-ray source. Comparison 
of the radio and optical properties of this extended radio emission argue against its association with 
either an H II region or supernova remnant. This is additional evidence that this ULX is not powered 
by a stellar-mass object whose emission is relativistically beamed towards the observer, and thus is 
either a super-Eddington source or intermediate mass black hole as suggested by optical observations. 
Implications of this result to future and existing radio studies of ULXs are discussed. 

Subject headings: galaxies: individual (Holmberg II) - stars: individual (Holmberg II X-l) - X-rays: 
galaxies - radio continuum: galaxies 


1. INTRODUCTION 

One of the more interesting discoveries of X-ray ob¬ 
servatories is that of unusually strong, compact non¬ 
nuclear sources in galaxies. These non-nuclear sources 
may have luminosities exceeding 2 x 10 39 ergs s _1 , mean¬ 
ing that should they result from isotropic emission at 
the Eddington limit the implied masses are larger than 
those expected for end products of normal ste llar evo¬ 
lution (about 20 M 0 ; iFrver fc Kalogeral 120011) . They 
might then be “intermediat e mass black holes” (IMBHs, 
IColbert fc Mushotzkvil 199?t) . an attractive possibility as 
these could fill in the black hole mass function between 
stellar-mass black holes and the supermassive (10 6 — 10 9 
M a ) black holes found in the centers of galaxies. Sim¬ 
ilarly, there is a large gap in the distribution of X-ray 
luminosities for known discrete sources between X-ray 
binaries (< 10 38 ergs s _1 ) and AGN (> 10 42 ergs s _1 ). 
Thus, ULXs are also often called intermediate X-ray ob¬ 
jects or IXOs. It is clear that if ULXs proved to be 
associated with IMBHs in the 10 1 2 — 10 4 M 0 range, they 
would no longer violate the Eddington limit. 

However, the identification of ULXs as IMBHs is only 
one possible theory. One obvious possibility is that ULXs 
are not physically associated with the host galaxies, be¬ 
ing background AGN (or foreground stars) whose lumi¬ 
nosity is thereby mistaken. While this may certainly 
explain some ULXs, the number of possible ULXs is well 
in excess of wh at wo uld b e exp ected based on random 
superpositions (iFahhiaunl ll 98911 . Supernova remnants 
(SNR) are capable of producing X-ray luminosities of 
~ 10 4 ° erg s s - 1 , and are viable exp lanations for some 
ULXs (e.g.- lFabian fc Terlevichll99fill . However, SNR do 
not exhibit the X-ray variability observed for most ULXs, 
which argues that ULXs are compact sources and thereby 
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presumably XRBs. 5 Thus, an attractive alternative are 
the more energetic hyp ernovae (potentially as sociated 
with gamma-ray bursts: fFrver X Wooslevlll 998 1. which 
could form both a com pact X-ray o bject and a large as¬ 
sociated optical nebula llWangi20f)21 . In fact, some ULXs 
have been demonst rated to be associated with large op- 
tical nebulae fe.g.. lWangll2004 iPakull X Mirion riQQi 


: (e.g IWanell2002t IPakull fc 

n i2nn.1t IKaaret et abll200l lU 


Roberts et al.ll2nn3t IKaaret et all 12004 iLehmann et alJ 

UmIT In many cases, the X-ray spectra of ULXs are 


consistent with expectati ons for accretin g black holes 
(e.g., rotating Kerr holes. (Makishima et a,l.ll20~00t) . How¬ 
ever, it is not required that these XRBs include a com¬ 
pact object more massive than a stellar mass black hole. 
iKing et all (l2nnifl proposed anisotropic emission from 
XRBs, thereby alleviating the need for super-Eddington 
emission. A similar model invokes relativistic beaming 
llKording et alJl2003l . with ULXs akin to microquasars 
whose jets are observed directly along our line of sight. 
There is also a model which allows super-Eddington 
emission fr om^th in accretion disks around stellar-mass 

We are engaged in a study of archival NRAO Very 
Large Array (VLA ) data which aims to identify ULXs 
with radio counterparts. Radio emission is an excellent 
tool for understanding ULXs and differentiating among 
the various progenitor models. Interferometric radio data 
provide good positional accuracy, uncomplicated by the 
dust extinction which often plagues optical analyses. Un¬ 
der the assumption that the radio emission is powered 
by the same source as the X-ray emission, the radio mor¬ 
phology can determine whether the emission is beamed. 
Unresolved, compact radio emission is consistent with 
beaming towards the observer where as extended emissio n 
rules out this possibility (e.g., as in IKaaret et alJ12003 ). 
Spectral indices garnered from multi-frequency observa¬ 
tions may also be assessed relative to known sources 
and emissio n mech ani sms (i.e., thermal emission or syn¬ 
chrotron) . iNeff et ab ({200581 noted the wealth of baseline 
expectations for the ratio of radio (at ~ 5 GHz) to X-ray 


5 Observed X-ray variability also largely precludes LILXs rep¬ 
resenting multip le u nresol ved XRBs, each of normal luminosity 
URoberts fc WarwickH lj^OOdl . 
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flux, enabling empirical comparisons between ULXs and 
more understood sources such as X-ray transients, SNR, 
star clusters, and AGN. Should some ULXs be associ¬ 
ated with IMBHs, this last comparison is particularly 
apt. Perhaps most importantly, radio data enable the 
computation of source total energies and lifetimes. These 
may be assessed in relation to specific ULX progenitor 
models as well as the evolution of such sources. 

Among the better studied ULXs is source X- 
1 in the dwarf irregular galaxy Holmberg II (or 
UGC 4305, DDO 50, and hereafter HoII) This ULX 
was f irst identified using the ROSAT HRI (|Zeza,s et all 
mm, and h as subsequently b een observed using 
both Chandra llK aaret et alJ 120041) and XMM-Newton 
I Dewangan et alJ I2004H . The proximity of HoII (3.05 
Mpc based on Cepheids, 3.39 M pc from HST obser¬ 
vations of the red gian t branch; iHoessel et alJ 119981 
iKarachentsev et al.ll2002l) has also made it a prime tar- 
get for optical observation. iPakull fc Mirionil (I2002T) 
noted He n emiss ion in the location of the ULX, 
which IK aaret et all 120041) confirmed using HST and 
iLehmamTet^iLr (l2004fl verified using several ground- 
based spectrographs. Both the X-ray data (in the spec¬ 
tral components and observed variation thereof) and the 
He II data argue that the ULX in HoII is associated with 
a black hole with L x > 4 x 10 39 ergs s -1 , and there¬ 
fore either a super-Eddington sour c e or an IMBH of a t 
least 25M 0 dDewangan et al.ll2004 IKaaret et al .1120041) . 
In the present paper, we discuss the radio properties of 
the ULX in HoII. 

2. DATA 

HoII has been the subject of several VLA observing 
campaigns. Our objectives are to detect radio emis¬ 
sion which is presumably faint, and do so at a resolu¬ 
tion high enough to resolve modest physical scales (1" = 
14.8 pc at the distance of 3.05 M pc assu med for HoII, 
chosen for consistency with IKaaret et ;)H 120041) . Thus, 
we sought longer duration observations performed in the 
larger VLA configurations. Two general programs fit 
these parameters, one at 1.4 GHz in the A configuration 
and the other at 4.86 GHz in the B configuration (each 
with a resolution of ~l , /5). There were also prior radio 
observations at ~15" resolut ion at 335 MHz, 1.4 GHz, 
and 4.86 GHz detailed in iTongue fc Westpfahll (Il995f) . 
These detected unresolved radio emission consistent with 
the ULX position, with corresponding flux densities of 
6.4±2.1 mJy, 1.37+0.08 mJy, and 0.77±0.05 mJy, re¬ 
spectively. 

The 1.4 GHz data were collected on 3 April 1994 as 
program code AT159, and consist of over 3 hours time 
on source. The raw data were obtained from the archive 
and loaded including antenna-specific weights. The flux 
scale was set using 3C286 with phase calibration provided 
by 0903+679 (J2000), following recommendations in the 
VLA Calibration Manual . 6 Imaging was done using the 
task IMAGR with multiple facets to alleviate problems 
such as the “3D effect.” The ( u,v ) data were gridded and 
tapered in a manner which resulted in a beam with good 
Gaussian characteristics, helping produce excellent sensi¬ 
tivity output maps. Multiple rounds of imaging and self 
calibration were performed until the final map was cre- 

6 \protect ittp://www.aoc.nrao.edu/~gtaylor/calib.html 
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Fig. 1.- 1.4 GHz map, with contours plotted at -2, 2, 3, 5, 8, 
and 13 times a base level of 17 pjy beam -1 . The restoring beam 
is shown at lower left. The Chandra ULX position is shown as a 
circle whose size represents its positional uncertainty. 


ated. This map has an rms noise of just 17 /iJy beam 1 
in the vicinity of the ULX position, where the beam size 
is T/9xl"5. 

Radio emission coincident with the ULX location is 
readily apparent (see Figure Q). Based on a Gaussian fit 
(using the task jmfit), its peak is located at 08 h 19 m 28!65 
+70°42'19" 0, about 1"6 f r om th e Chandra ULX position 
reported in IKaaret et al. l|2004lh The error in the radio 
position is less than O'/l, although the radio emission is 
clearly extended and covers the Chandra position. Based 
on the surface density of detected radio sources within 
the central 10' of the radio image, the probability of a 
chance ULX-radio superposition within 1"6 is well under 
0.1%. The radio emission has deconvolved major and 
minor axes of 3"7 (lcr range of 3"4 - 4"0) and 2"7 (2"4 
- 2'.'9) at a position angle of 91° (80° - 104°). The flux 
is 1.174+0.085 mJy, after correction for primary beam 
attenuation. Emission on slightly l arger scales is poss ible 
given the 1.4 GHz flux reported in lTongne Kr. West.nfahll 
f 1995 + although the two fluxes do not differ significantly. 

The 4.86 GHz data consist of three related programs 
totalling four observations. Unfortunately, each observa¬ 
tion is relatively short at ~T0 minutes time on source. 
Program AR165 observed HoII on 18 January 1988 and 
14 March 1989, while program AR227 and AR258 ob¬ 
served it on 29 July 1990 and 6 December 1991, respec¬ 
tively. Each observation was performed with effectively 
identical observational setups, making it fairly easy to 
combine the data into a single ~ 40 minute observation 
as well as analyzing each individually. 

Data reduction was achieved using the same basic pro¬ 
cedures as the 1.4 GHz data except in this case the phase 
calibrator was 0642+679 and the lack of strong sources 
within the fields prevented self-calibration. The ULX 
was not significantly detected in any of the individual 
observations, with rms noises in the ULX vicinity rang¬ 
ing from 58 to 79 p,Jy beam -1 . Peaks in the radio emis¬ 
sion within 2" of the ULX Chandra position were de¬ 
termined, and although the flux between the different 
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Fig. 2.— 4.86 GHz map. All contours and symbols the same as 
Figure IT1 with the exception of the base contour level being 44 /ijy 
beam - . 



circle sized to represent the positional uncertainty. The H II region 
extends well to the west of this figure. 


observations differed by a factor of two this variability 
is not statistically significant. At this point, the ( u,v ) 
data for the four observations were combined, weight¬ 
ing the individual observations by the inverse of their 
rms noise squared. The rms noise of the final map cre¬ 
ated using the combined data was 44 /rJy beam -1 with 
a beam size of l'/7xl'/5. Weak extended emission having 
an integral flux of 0.677±0.207 mJy covers the ULX po¬ 
sition (see Figure |3l. The deconvolved size is 374 (273 - 
475) by 179 (170 - 2") with a position angle of 37° (10° - 
66°). The peak in this emission, located at 08 h 19 m 28!84 
H-70°42 , 1975 (accurate to ~072 and within 077 of the 
Chandra position), does not constitute a traditional 5cr 
detection. However, imaging runs using pure natural 
weighting of the (u,v) data with a slight taper were able 
to yield detections at over 6er significance. The inte¬ 
gral flux of the source was consistent in all imaging runs 
and add itionally is consiste nt with the value reported in 
iTongue fc Westofahl (I1995T) . 

Taken together, the 1.4 GHz and 4.86 GHz data imply 
a spectr al in dex of a = 0.4 4 ± 0.31 (where S v oc u~ a ). 
ITongue fc Westofahl ( 1995 7) determined a spectral index 
of 0.5±0.1 over the same frequency range for the region 
coincident with the ULX, based on their lower resolution 
data. The greater uncertainty in our value is simply the 
result of the fairly noisy 4.86 GHz data. 

3. DISCUSSION 

The location of the ULX in HoII has been shown to lie 
within extended radio emission. Should the radio emis¬ 
sion arise from the same physical source as the X-ray 
emission, the resolved nature of the radio argues against 
relativistic beaming along the line of sight. In fact, the 
size of the radio emission implies a region of ~ 60pc by 
~40pc, elongated in the East-West direction. The east¬ 
ern region of the radio e mission is coi ncident with the 
He II nebulae reported in lKaaret, et ah l (|2004fh as shown 
in Figure 0 

The obvious question here is to what extent the ra¬ 
dio emission is linked to the ULX. The ULX position 


is associated with a region of star for mation indicate d 
by Ha and other emission lines ( e.g.,lPu che et al1ll992t 
iKaaret et alJ 12004 iLehmann et al .1 12004 . and H II re¬ 
gions are known sources of radio emission through ther¬ 
mal bremsstrahlung and non-thermal emission arising 
in supernovae. The derived radio spectral index is 
marginally steeper than that of bremsstrahlung which 
has a ~ 0.1, suggesting that all the emission is not sim¬ 
ply due to this. We can also get some idea of the contri¬ 
bution of bremsstrahlung to the net radio fl ux thro ugh 
t he r eported emission line measurements. IKaaret, et a.U 
d2QQ4|) find a net flux of 1.1 x 10~ 14 ergs s _1 for H/3, and 
ILehmann et all ( 2004 1 find that the H II region has little 
extinction (i.e., the ratio of Ha to H/3 is close to 2.85). 
For thermal b remsstrahlung emission in H II regions, 
iMoorwood fc Salinaril (1983) derive the expression: 

'*• - = 271 x 10 ^w rOJX( 5 

Using a line intensity ratio of 0.0834 for j j] Ta /j (e.g., 

Osterbrock for T = 10 4 K), this predicts a 4.86 GHz 
flux of only 0.034 mJy. Thus, the likely contribution 
of thermal bremsstrahlung to the total radio emission 
is quite small, about 5%. Qualitatively, we also note 
that the radio and He II emission appear only at the 
eastern edge of the full H II region, which is ~ 15" in size. 
Furthermore, the He n emission line is narrow and lacks 
the associated “bump” from nearby N III, C III, and C IV 
seen in Wolf-Rayet spectra. That the line is narrow also 
shows the He II is optically thin, further underscoring 
that the calculated lower bound on the X-r ay luminosity 
(4 x 10 39 ergs s _1 ) of the source reported in lKaaret et alJ 
( 2004 *1 is indeed a true lower limit. 

Alternatively, we may consider SNR. The total 
4.86 GHz flux is remarkably comparable to what would 
be emitted by Cas A were it to lie within HoII, and 
the derived spectral ind ex is consistent with s hell-type 
supernovae (a = 0.45. iClark fc Caswelllll976 *). How¬ 
ever, the radio morphology does not exhibit a shell- 
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type appearance but ra ther is filled, or “plerionic” (e.g., 
IWeiler & PanagiJll978 1. These types of SNR have flat¬ 
ter spectral indices (a ~ 0 ) and are of much smaller size. 
In fact, the radio-emitting region around the HoII ULX 
is very large for its lumino sity, placing it as an outlier 
to any X — D relation (e.g., Similarly, the 

integral field spectroscopy of lLehmann et all ( 1200 41 in¬ 
dicates that the line emission is not that expected for 
supernovae nor i s the X-ray emission co nsistent with an 
SNR origin (e.g., iDewangan et al.ll2004l) . 

The radio data also allow for calculation of total energy 
and cooling time un der the assumption of equipartition 
(e.g., lLongaiilll994 j. For this calculation, we made the 
simple approximation that the radio source was a sphere 
whose projected area is equivalent to the measured size 
at 1.4 GHz, and that the filling factor within this sphere 
was unity (the effect of lowering the filling factor is to 
reduce the minimum energy and cooling time). The to¬ 
tal radio luminosity of the source was estimated by as¬ 
suming a power law between the observed frequencies 
of 1.4 GHz and 4.86 GHz. The resulting minimum en¬ 
ergy was 2.6 x 10 49 77 4 / 7 ergs, implying a cooling time of 
2.5 x 10V / 7 years (where 77 is the ratio of baryon to 
electron particle energy). The calculated minimum mag¬ 
netic field strength was 13?7 2 / 7 /iG. Derived cooling times 
under the assumption that the radio emission is thermal 
bremsstrahlung are about an order of magnitude larger. 

The collection of X-ray, optical, and radio data show 
that the radio emission is not due to a simple SNR or 
H 11 region, and that the radio emission is positionally 
associated with the X-ray source and the He 11 nebula. 
We are thus led to the conclusion that the radio emis¬ 
sion and the ULX are causally connected. In searching 
for other examples of accreting binaries with extended lu¬ 
minous radio emission, we note some similarities to the 
Galactic system of the W50 nebula and the jet source 
SS 433. In each case, the radio emission covers an extent 
of ~ 50 pc with a t otal spectral index of about 0.5 (e.g., 
iDubner et all 11998 1. However, the HoII nebula is over 
an order of magnitude more luminous and correspond¬ 
ingly implies a much greater energy in relativistic parti¬ 
cles (5 x 10 4 y/ 7 ergs for W50-SS 433 vs. 2.6 x 10 49 ?7 4 / 7 
ergs for HoII: IDubner et aHll998 ). 

Computation of a cooling time also enables an in¬ 
vestigation o f poss ible source motion as discussed in 
iKaaret, et a.1.1 (2 004 1. The simple hypothesis is that ULXs 
are associated with stellar clusters, and their X-ray emis¬ 
sion declines as they move away from these clusters. The 
HoII ULX is bright, but not close to any apparent stellar 
cluster. Thus, for this model to be correct the HoII ULX 
would need to be travelling at high velocity. Using the 
major axis size of ~ 60 pc and the cooling time of 2.5 x 10' 
years, the velocity is only about 2 km s -1 . This is much 
less than the expected kick velocities of X-ray binaries 
(unless there is a very substantial radial component to 
the velocity), and would seem to preclude the ULX be¬ 
ing a young system ejected from a stellar cluster at high 
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velocity. 

What are the prospects for detecting similar ULXs 
in t he radio? B ased o n the X-ray luminosites reported 
in IKaaret et alJ [ 2004 1 and our radio results, a repre¬ 
sentative value of R\ (the radio to X-r ay ratio, see 
iTerashima fc Wilsonll2003t iNeff et, aHl2003H for the HoII 
ULX is 5 x 10~ 6 (about the same as for the NGC 5408 
ULX IKaaret et all2003tfNeff et alJ2003h . Assuming this 
Rx and Lx = 10 39 ergs s -1 , a one-hour integration with 
the VLA can detect such a source at the 5er level provided 
it lies within a distance of 2.8 Mpc and is unresolved. If 
this hypothetical source has an intrinsic size of 40 pc 
(hence about 3" across), the one-hour observation will 
fail to detect it unless the VLA is in its more compact C 
or D configurations. 

The above explanation fits t he re s ults of two other ra¬ 
dio studies of ULXs. IKaaret et al. (l2003t) used archival 
ATCA data to detect the ULX in NGC 5408 at 4.8 GHz 
(where it was unresolved) but not at 8.64 GHz. The 
resolu tions were 2 "9 x 1"8 and 1"6 x 1"0, respectively 
(see iStevens et al.ll2002lh At the assumed distance of 4.8 
Mpc, 1" = 23 pc. Given the noted 40 /iJy beam -1 rms 
of the 8.64 GHz data, emission extended over ~40 pc 
might easily be missed. This scenario is roughly consis¬ 
tent with the spectral index of the HoII data, provided 
the 4.8 GHz NGC 5408 emission is only marginally un¬ 
resolved. This theory will be testable in the near future, 
as an NGC 5408 monitoring campaign using the VLA is 
currently in pro gress (Kaaret, private communication). 
Similarly. iKording et al. ( l2005t) performed a monitoring 
campaign of nine nearby ULXs (within 5.5 Mpc) with 
possible detections of two compact sources, believed to 
be SNR. The observations were performed with the VLA 
primarily in its largest configuration (A) at 8.46 GHz, 
with a resolution of about 0"25. Their excellent sensi¬ 
tivities of ~ 15 /rJy beam -1 are still nearly an order of 
magnitude larger than necessary to detect any diffuse 
emission akin to that observed for HoII. 

In summary, we have identified an extended radio 
source coincident with the ULX in HoII. The peak of the 
radio emission lies within l'/6 of the ULX location, im¬ 
plying a very small probability that the radio and X-ray 
sources represent a random superposition. The proper¬ 
ties of the radio emission argue against simple interpre¬ 
tations as an H II region or SNR, particularly when con¬ 
sidered in light of the He n nebula and attendant optical 
spectroscopy. We believe that this association, very dif¬ 
ferent from what has been observed for accreting X-ray 
sources within the Galaxy, furthers the identification of 
HoII X-l as an intermediate mass black hole. 


The authors thank Michael Rupen for supplying the 
4.86 GHz radio data, Schuyler Van Dyk for a preliminary 
1.4 GHz image, Philip Kaaret for the He II image, and 
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